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Abstract 

We investigate the branching structure coded by the excursion above 
zero of a spectrally positive Levy process. The main idea is to identify the 
level of the Levy excursion as the time and count the number of jumps 
upcrossing the level. By regarding the size of a jump as the birth site of 
a particle, we construct a branching particle system in which the particles 
undergo nonlocal branchings and deterministic spatial motions to the left 
on the positive half line. A particle is removed from the system as soon as it 
reaches the origin. Then a measure-valued Borel right Markov process can 
be defined as the counting measures of the particle system. Its total mass 
evolves according to a Crump-Mode- Jagers branching process and its sup- 
port represents the residual life times of those existing particles. A similar 
result for spectrally negative Levy process is established by a time reversal 
approach. Properties of the measure-valued processes can be studied via 
the excursions for the corresponding Levy processes. 

AMS 2010 subject classifications: Primary 60J80, 60G51; Secondary 60J68. 

Keywords: Levy process, spectrally one-sided, subordinator, branching par- 
ticle system, non-local branching, Crump-Mode- Jagers branching process. 

1 Introduction 

Branching processes embedded in processes with independent increments have been stud- 
ied by many authors. The study yields detailed information and understandings in the 
two classes of processes. In particular, Dwass [7j constructed branching processes from 
simple random walks. To study random walks in random environment Kesten et al [H] 
constructed a Galton- Watson process with geometric offspring law from a simple random 
walk. Multitype branching processes have also been introduced in the study of random 
walks in random environment; see JEHHI1II2] and the references therein. Since continuous 
state branching processes and Brownian motions arise as the scaling limits of Galton- 
Watson processes and simple random walks, respectively, we may naturally expect some 
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branching structures embedded in a Brownian motion. The well-known Knight-Ray the- 
orem brings an answer to this question; see also [16] [21] . 

Le Gall and Le Jan jTH [19] recovered a deep connection between general continuous 
state branching processes and spectrally positive Levy processes. Furthermore, Duquesne 
and Le Gall [4, 5J showed that the branching points of a Levy tree constructed in [TS] are 
of two types: binary nodes (i.e. vertex of degree three), which are given by the Brownian 
part of the Levy process, and infinite nodes (i.e. vertex of infinite degrees), which are 
given by the jumps of the Levy process. The size of the jump is also called the size of the 
corresponding infinite node (or the mass of the forest attached to the node). 

In the interesting recent work [14], Lambert used spectrally positive Levy processes 
for the first time to code random splitting trees. In the population dynamics represented 
by the splitting tree, the number of individuals evolves according to a binary Crump- 
Mode- Jagers process. It was proved in |T3] that the contour process of the splitting 
tree truncated up to a certain level is a spectrally positive Levy process reflected below 
this level and killed upon hitting zero. From this result Lambert derived a number of 
properties of the splitting tree and the Crump-Mode- Jagers process. 

The purpose of this paper is to give a formulation of the branching structures of 
spectrally one-sided Levy processes in terms of measure-valued processes, which we call 
single-birth branching particle systems. Those structures are undoubtedly conveyed by 
the random splitting trees, so we could have derived the results from those of Lambert 
[T4] . However, we think a simple construction of the branching particle systems directly 
from the Levy process is of interest. In addition, we show that the branching systems are 
Borel right Markov processes in a suitable state space and characterize their transition 
semigroups using some simple quasi-linear integral equations. Those properties make 
the branching systems easier to handle than the Crump-Mode- Jagers processes. A more 
precise description of the branching structures is given in the next paragraph. 

Let us consider a typical trajectory of the spectrally positive Levy process with nega- 
tive drift {S t : t > 0} started from a > and killed upon hitting zero; see Figure 1'. Let 
{yi : i = 1, 2, 3} denote the sizes of jumps. Then the sample path of a branching particle 
system can be obtained in the following way: At time zero, an ancestor starts off from 
a > and moves toward the left at the unit speed. At times z\ and z 3 , it gives birth to 
two children at positions y\ and y%, respectively. At time Z2, the first child of the ancestor 
gives birth to a child at position y2- Once an individual hits zero, it is removed from the 
system. So the ancestor dies at time a and its two children die at times zi + yi and z 3 + y 3l 
respectively. 

From the structures described above, we use a time reversal to derive a similar result 
for spectrally negative Levy processes with positive drift. We will see that the branching 
systems we encounter here are actually very special cases of the models studied in [3j [20] . 
Unfortunately, by now we can only treat Levy processes with bounded variations as in [H] . 
An interesting open question is to give a description of the branching structures of general 
spectrally one-sided Levy processes in terms of measure-valued branching processes. We 
hope to see the precise formulation of such structures in the future. 
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time of branching process 
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Figure 1 

The rest of this paper is arranged as follows. In Section 2, we introduce some branching 
particle systems on the positive half line involving nonlocal branching structures. In 
Section 3, we extend the model to the case with infinite branching rates. In Section 4 the 
result on the branching structures in spectrally positive Levy processes with negative drift 
is established. In Section 5, we derive the branching structures for spectrally negative Levy 
processes with positive drift by a time reversal approach. Some properties of our branching 
systems are studied in Section 6. In Section 7, we discusses briefly the connection of the 
branching systems with the Crump-Mode- Jagers models. 

Notations. Write R + = [0, oo). Given a metric space E, we denote by B(E) the 
Banach space of bounded Borel functions on E endowed with the supremum/uniform 
norm "|| • ||". Let C(E) be the subspace of B(E) consisting of bounded continuous 
functions on E. We use the superscript "+" to denote the subset of positive elements of 
the function spaces, e.g., _B + (1R + ) and C + (0,oo). Let M(E) denote the space of finite 
Borel measures on E endowed with the topology of weak convergence. Let N(E) be the 
set of integer- valued measures in M(E). For a measure fi and a function / on E write 
(//, f) = J fdfi if the integral exists. Other notations will be explained when they first 
appear. 

2 Branching systems on the positive half line 

We begin with the description of a branching system of particles on IR + . Suppose that 
a > is a constant, r\ = rj(dx) is a probability measure on (0, oo) and g = g(z) is a 
probability generating function with g'(l) < oo. Let {£ t : t > 0} be the M + -valued 
Markov process defined by £ f := (£ — £) V 0. Let F(Q, •) be the unit mass at S G iV(R + ). 
For x > 0, let F(x, •) be the distribution on iV(R+) of the random measure 

z 

i=l 

where Z is an integer-valued random variable with distribution determined by g = g(z) 
and {Yi,Y 2 ,- ■ ■ } are i.i.d random variables on (0, oo) with distribution rj(dx). Here we 



assumed that Z and {Y\, Y 2 , ■ ■ ■ } are independent. 

Suppose that we have a set of particles on K + moving independently according to the 
law of {£t : t > 0}. A particle is frozen as soon as it reaches zero. Before that at each 
a-exponentially distributed random time, the particle gives birth to a random number of 
offspring according to the law specified by the generating function g = g(z), and those 
offspring are scattered over IR + independently according to the distribution r](dx). It 
is assumed as usual that the reproduction of different particles are independent of each 
other. Let X t (B) denote the number of particles in the set B G ^(R + ) at time t > 0. By 
Dawson et al [3, p. 103] one can see that {X t : t > 0} is a Markov process on iV(IR + ) with 
transition semigroup (Qt)t>o defined by 



>,/> 



Q t (ji,du) = exp{-{fi,U t f)}, f G B + (R + ) 



(2.1) 



AT(0,oo) 



where (t,x) h-» U t f(x) is the unique positive solution of 



-Utf(x) 



e -/(0»-t)vo) _ a 



,-u t -.f«x-.)va) d8 + a / e -K-J(o) 1{i<s}(is 
Jo 

+ a I e- Dt -^ x - s h {x>s} g(( V ,e- Dt -° f ))ds; 



see also Dawson et al [3l pp. 95-96] and Li [201 P-98]. By Proposition 2.9 of [20], the above 
equation can be rewritten as 



e -Utf(x) _ e -/((*-t)V0) _ a 



)l { , >a} [l -^, e -^-/»]d s . (2.2) 



By Proposition A. 49 of [20], for / G 5(R + ) there is a unique locally bounded solution 
(t,x) h-> 7Ttf{ x ) to the equation 



(2.3) 



7f t /(x)=/((x-t)V0)+a^(l) / l {!e >.}fa,7ft-./>ds. 

Moreover, the linear operators (vf*)t>o on 5(R + ) form a semigroup and 

IMI < \\f\\e a9 ' (1)t , t>0. 



(2.4) 



Proposition 2.1 For t > and / G 5 + (lR + ) we nave U t f < n t f and for t > and 

/ G 5+(M + ) we have 



(vJ)Qt{^dv) = (n,ir t f)- 



-) 



(2.5) 



Proo/. For t > and / G B + (M+) one can use and (Q to see 
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6»=0 



Then (2.5) follows by differentiating both sides of (2.2). By (2.1), (2.5) and Jensen's 
inequality it is clear that U t f(x) < Tt t f{x) for x > 0. By linearity we also have (2.3) and 
dZ5b for / G B(R+). □ 
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Proposition 2.2 For any f 6 themappingt i-)- U t f(-+t) from [0, oo) toi? + (]R + ) 

is increasing and locally Lipschitz in the supremum norm. Moreover, for any t > r > 
we have 



< e -&r/(*+r) _ e -U t f(x+t) < a(t _ r) _ 



(2.6) 



Proof. For any > one can use (2.2) to see 



e -u t f( x +t) = e -f(xvo) 



ale 
'o 



'^ +s )l { , +s> o } [l-^((r/,e-^))]^. 



Then we have (2.6). Since t i— )■ i7 t / is locally bounded by Proposition 



2.1 



£/*/(• + i) is increasing and locally Lipschitz in the supremum norm. 



we see t \-> 

□ 



Proposition 2.3 For any f e 5 + (M + ) the function (t,x) h-> U t f(x) is the unique locally 
bounded positive solution of 

UJ(x) = f((x-t)V0) + a [ 1 {X> , } [1 -g((n,e-°<->f)j\ds. (2.7) 



Proof. For notational convenience, in this proof we set f(x) = /(0) and Utf(x) = Utf(0) 
for all x < and t > 0. Let = t < t\ < ■ ■ ■ < t n = t be a partition of [0, t]. For x G M + , 
we can write 



U t f(x) = f\ x -t) + Y J \Ut-u_J\x - U. x ) - U t - ti f(x - U) 



{21 



i=l 



Note that Proposition 2.3 implies U t „ ti _J(x - t^) - U t „ t J{x - U) > 0. By (2.2), (2.6) 



and Taylor's formula, as t$ — £j_i — >■ 0, 



Ut-ti-Jix-U^-Ut-ttfix-ti) 



3 ft-t i _ 1 /(a'-ti-l) 



1 



-U t -tJ(x-h) _ -Ut-t^Jix-ti-!) 



[l + £i(s, x)] l {x - u _ 1>s } [l - #((?7 



g-ft-t i _ 1 -« 



+ o(tj - ti_l) 

/ ))]ds + o(t i -t,_ 1 ), 



where 



£l -(s,x) = e ^-*i-i/(»-*i-0 



-Ct-t^i-s/C^-ti-i-s) _ p-Ut-ti^fix-ti-i) 



By Propositions |2.1| and |2.2| one can see that 

0<ei(s,z) < a^-ti-Oexpdl/He^W*}, < s^^Vi. 
It then follows that 

Ut-u.Jix - t 4 _i) - U t -tJ(x - U) 



ti—ti-1 



Ux-t i - 1 >s} [l-g((v, 



-Ut- 



* f ))]ds + o(ti-t 



i-l, 



Mx>s } [1 - g((v, e- 0t -° f ))]ds + o{U - U-x). 



Substituting this into (2.8) and letting maxi<j<„(tj — ti-i) — > we obtain <\2.7 ). The 
uniqueness of the solution of the equation follows from Proposition 2.18 in (20]. □ 



/ e~^Q t {^dv) = e-^\ 

JN(0,oo) 



Theorem 2.4 There is a Borel right transition semigroup (Qt)t>o on N(Q, oo) defined by 

feB+(0,oo), (2.9) 
where (t,x) i-» U t f(x) is the unique locally bounded positive solution of 

U t f(x) = f(x-t)l {x>t} + a f l {x> ^ s} [l-g((r),e- u *f))]ds, t>0,x>0. (2.10) 



Proof. It is not hard to see that (2.10) is a special cases of (2.21) in j2(Jl p. 39]. By (2.2) 
we have Utf(0) = /(0) for all t > 0. Conseq uentl y, if { X t : t > 0} is a Markov process 
with transition semigroup (Qt)t>o defined by (2.1) and (2.7), then {X t \(o t00 ) : t > 0} is a 



Markov process in N(0, oo) with transition semigroup (Qt)t>o defined by (2.9) and (2.10). 



By Theorem 5.12 of [20], we can extend (Qt)t>o to a Borel right semigroup on the space 
of finite measures on (0, oo). Then (Qt)t>o itself is a Borel right semigroup. □ 



By Proposition 2.1 we have the following: 



Proposition 2.5 For every f G 5(0, oo) there is a unique locally bounded solution 
(t,x) h-> ir t f(x) of 



(2.11) 



7T t f(x) = f((x - t) V 0) + ag'(l) / l {x>s} (v, n t _J)ds. 

Jo 

Moreover, the linear operators (irt)t>o on B(0, oo) form a semigroup and 



{u, f)Q t (fi, du) = (fi, 7T t f), t>0,fe B(0, oo) 



(2.12) 



JV(0,oo) 



Proposition 2.6 We have U t f(x) < n t f(x) < ||/||e a9 ' (1)i for t > 0, x > and f e 
S(0,oo). 



A Markov process in A^(0,oo) with transition semigroup (Qt)t>o defined by (2.9) and 



(2.10) will be referred to as a branching system of particles with parameters (g,a,rj), 
where g is the generating function, a is the branching rate and r] is the offspring position 
law. 



3 The system with infinite branching rate 

In this section, we consider a system of particles, which can be thought of as a branching 
system with infinite branching rate. Let p(x) = x for x G (0, oo). Let B p (0, oo) be the set 
of Borel functions on (0, oo) bounded by p ■ const. Let C p (0, oo) be the subset of B p (0, oo) 
consisting of continuous functions. Let M p (0, oo) be the set of Borel measures p on (0, oo) 
satisfying (p,p) < oo. Let N P (Q, oo) be the set of integer-valued measures in M p (0, oo). 
We endow M p (0, oo) and N p (0, oo) with the topologies defined by the convention that 

fjL n ->■ fi if and only if (p n , f) -)> (p, f) for all / G C p (0, oo). (3.1) 

We say a function (t,x) !->■ Wt(x) on [0, oo) x (0, oo) is locally p-bounded if 

sup sup \p(x)~ u s (x)\ < oo, t > 0. 

0<s<t xe(0,oo) 

Let c > be a constant and let Il(<iz) be a er-finite measure on (0, oo) such that 
(II, p) < c. Given / G B+(Q, oo), we consider the following evolution equation: 



Utf(x) = f(x - t)l {x>t} + c 1 / l {;r>s }(n, 1 - e 

J o 



)ds. (3.2) 



Lemma 3.1 For each f G 5^(0,00) tiere is at most one locally p-bounded positive 
solution of (3.2). 

Proof. Suppose that (t,x) y-¥ U t f(x) and (t,x) V t f(x) are two locally p-bounded 
solutions of (3.2). Let 

l T (x)= sup \p(x)- 1 (U t f(x)-V t f(x))\. 

0<t<T 

Then for any < t < T we have 



\U t f(x) - V*/(ar)| < c" 1 / l {a;>s} (n, le"^"^ - e' v ^\)ds 



< c- 1 f 1 {X>S} (TI, |£/ t _J - 7 t _ a /|)ds 

< c- 1 / l w d S ||Z T ||(n,p) < c-V(^)||/r||(n,p), 

Jo 



which implies \\It\\ < c 1 ||/t|| (n, p). Then we have ||/(T)|| = as (II, p) < c. □ 



Proposition 3.2 For eacn / G l?j~(0,oo), there is a unique locally p-bounded positive 



solution (t,x) i — y U t f(x) of (3.2) and the solution is increasing in (II,/) G M p (0, oo) 



x 



B+(0, oo). Furthermore, the operators (U t )t>o on _B+(0, oo) form a semigroup and 



Ip-'ujw^ic-i^p^Wp-'fi 



t > 0. 



(3.3) 



Proof. Step 1) We first assume that II G M(0, oo) and / G -B + (0, oo). By Theorem 4.1 
there is a unique locally bounded positive solution (t,x) U t f(x) of (3.2). This solution 
can also be constructed by a simple iteration procedure. In fact, if we let Uo(t, x) — and 
define u n (t,x) = u n (t,x, f) inductively by 

pt POO 

u n {t, x) := f(x - t)l [x>t} + c" 1 / l {x>s} ds / [1 - e-^-^-^Uidz), (3.4) 

JO JO 

then u n (t,x) — )■ Utf(x) increasingly as n — >• oo; see Proposition 2.18 of [20J. Using this 
construction one can see that the solution of (3.2) is increasing in (II,/) G M(0, oo) x 
5+(0,oo). 

Step 2) Next, we assume that II G M(0, oo) and / G -B+(0, oo). Let f k — f A k for 
> 1. Let (i, x) h-> U t fk(x) be the unique locally bounded positive solution of (3.2) with 
/ replaced by fk- According to the argument above the sequence {U t fk} is increasing in 
k > 1. By (3.2) and Proposition 2.6 we have 



Utfk(x) < \\p- L f k \\p(x) + c 



-i 



/ £/ t _ 8 / fc (z)n(dz) 



< 



/o- i AII + c- 1 ||MI(n,i)exp{ c - i (n,i)t} P (x 



Thus (t,x) i — y U t fk{x) is locally p-bounded. On the other hand, if we set 



l k {t,x) := sup U s f k (x) 

0<s<t 



then 



h(t,x) < \\p 1 f k \\p(x) + c l p{x) sup 

0<s<t J(0,oo) 

|p- 1 /||+c- 1 ||p- 1 /*(t)||<n ) p>lp( a :). 



U s f k {z)U{dz) 



< 



It follows that 



which implies 



p(x)-%(t,x) < Wp-ifW+c-^p-XWUlLtp), 



\\p-%(t)\\ < 



p VII 



c||p-7II 



(3.5) 



i-c-i(n, P ) c-(n, P )' 

In particular, we have 

\\p- l U t f k \\ <c\\p- l f\\{c-{Xl,p))-\ t>0. 

Then the limit U t f(x) := lim^oo U t f k {x) exists. It is easy to see that (t,x) i— )■ U t f(x) is 
a locally p-bounded positive solution of (3.2) satisfying (3.3). 

.Step 3) In the general case, let Hk(dz) = l{ z >i /k}H(dz) for k > 1. For / G -B + (0, oo) 
let (t, x) I—?- 

C/ t (fe) /(^) be the unique locally p-bounded positive solution of (|3.2[) with 
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II replaced by IT^. By the second step, we can define U[ f by the equation for any 
/ G Bp(0, oo). The sequence {U^ f} is increasing by the first and the second steps. 
As in the second step one can see the limit U t f(x) := lim^oo fix) exists and is a 



locally p-bounded positive solution of (3.2) satisfying (3.3). The uniqueness of the solution 
follows from Lemma 3.1, which yields the semigroup property of {U t )t>o- D 



Proposition 3.3 For each f G B p (0, oo), there is a unique locally p-bounded solution 

(t,x) I y 7T t /(x) Of 



(3.6) 



■Ktf(x) = f(x-t)l {x>t} + c / l {x>t _ s} (TL, 7r s f)ds 



Furthermore, the solution is increasing in (II, /) G M p (0, oo) x B p (0, oo) and (7T t ) t > is a 
semigroup of linear operators on B p (0, oo) such that 



p-^ t f\\<{c-(Tl,p))^\\p-^l 



t > 0. 



(3.7) 



Proof. For / G B+(0, oo) one can obtain (3.6) by differentiating both sides of (3.2), and 



(3.7) follows by (3.3). By the linearity, the equation has a solution for any / G B p (0, oo) 
and (3.7) remains true. By Proposition 3.2 one can see the solution is increasing in 



(LT, /) G M p (0, oo) x B p (0, oo). The uniqueness of the solution follows by a modification 
of the proof of Lemma |3.1[ □ 



Theorem 3.4 There is a Borel right semigroup (Qt)t>o on N p (0, oo) defined by 

e-MQ t (n,dv) = e-W*f) f f G R+( , oo), (3.8) 



/ 



AT p (0,oo) 



where (t, x) >-)■ U t f(x) is the unique locally p-bounded positive solution of (3.2). Further- 
more, we have 



(f, f)Qt{v, dv) = (/!, n t f), f G B p (0, oo) 

'Af p (0,oo) 

where (t, x) i— )■ ittf{%) JS the unique locally p-bounded solution of 



(3.9) 



TTt/C^) = /(^ - *)l{a:>t} + C 1 / l{ x>t - s }(U, 7T s f)ds. 

JO 



(3.10) 



Proof. Let 



Theorem 



2.1 



i>o be defined as in the last step of the proof of Proposition 3.2 By 
we can define a Borel right semigroup (Qt )t>o on N(0, oo) by 



-<M,C/ t (fc) /> 



/G5 + (0,oo). 



(3.11) 



AT(0,oo) 



In view of (|2.12|) and (|3.7|), if p G N p (0,oo) is a finite measure, we can regard Q[ k \p 



as a probability measure on N p (Q, oo). Clearly, N p (0, oo) is a closed subset of M p (0, oo) 
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and the latter is an isomorphism of M(0, oo) under the mapping v(dx) t— > xvidx). By 
Theorem 1.20 of [20J and the last step of the proof of Proposition 3.2 one can see ( 3.8[ ) really 
defines a probability measure ■) on iV p (0, oo) for any finite measure fi 6 iV p (0, oo). 

By approximating fi G iV p (0, oo) with an increasing sequence of finite measures, we infer 



the formula defines a probability kernel on iV p (0, oo). Here (3.2) can be regarded as a 
special form of (6.11) in [20]. By Theorem 6.3 in [20J, we can extend (Qt)t>o to a Borel 
right semigroup on M p (0, oo). Then we infer that (Qt)t>o itself is a Borel right semigroup. 



The moment formula ( |3.9 ) can be obtained as in the proof of Proposition 2.1 



□ 



A Markov process in N p (0, oo) with transition semigroup (Qt)t>o defined by (3.2) and 



(3.8) will be referred to as a single-birth branching system of particles with offspring 
Clearly, when II is a finite measure on (0, oo) 



position law II 
case of the model introduced in the last section. 



this reduces to a special 



4 Subordinators with negative drift 

In this section, we give a description of the branching structures in subordinators with 
negative drift. Set 

C 1 (M) = {/ G C(M) : f is different iable and has bounded derivative.} 

Let c > be a constant and let IT be a cr-finite measure on (0, oo) satisfying (II, p) < c. 
Suppose that {S t : t > 0} is a subordinator with negative drift generated by the operator 
A given by 

POO 

Af(x)= [f(x + z)-f(x)]IL(dz)-cf(x), feC'iR). (4.1) 

JO 

We assume S = a > 0. Our assumption implies that S t — > — oo as t — > oo, so the hitting 
time 

To := inf{t > : S t < 0} 

is a.s. finite. For t > set 

J(t) := {u e [0, r "] : S u _ < t < S u } (4.2) 
with the convention that £o_ = 0. Then we define the measure- valued process 

X t = 5s u -t, t > 0. (4.3) 
ueJ(t) 

It is easy to see that X = 5 a . 



Theorem 4.1 The process {X t : t > 0} is a single-birth branching system in N p (0, oo) 
with transition semigroup (Qt)t>o defined by (3.2) and (3.8). 
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Proof. Step 1 ) We first assume U(dz) is a finite measure on (0, oo). In this case we clearly 
have P{#J(i) < oo for all t > 0} = I. Let 

C(t) = {u G [0, Tq] : S u = S u ^ = t} and ((f) = #C(t). 

We can write C(t) = {ti(£), • • • , T((t)(£)} by ranking the elements in increasing order. Let 
T (t) = and 

<7i(t) = inf{n > ri_i(t) : S„ > £}, z = 1, 2, - • • , C(t). 
Then it is easy to see that J(t) = {ai(t), ■ ■ ■ , cat)(t)} and 

i=i 

In particular, we have Q(t) = #J(t). Write M t = min < r <t S r and L t = S t — M t . Set 
1] Q = and for > 1 define inductively 

C fc = inf{t > : $ ^ M t } and % = inf{t > ( k : S t = M t }. 



time of X 




time of 5 



Cs % Ca 

Figure 2 



It is clear that a — S^- is an exponentially distributed random variable with parameter 



A/c, where A 
e k (t) = (L Ck+t 



11(0, 00). By the memoryless property one can see S, 



Vk-l 



- S t 



5*^_ is also exponentially distributed with parameter A/c for each k > 1. 



Let 



-L^ fe „)l{ t<r?fc _^ fc } and let F(dw) denote the distribution of {e\{t) : t > 0} on 



_D + [0, 00), the space of positive cadlag functions on [0, 00) equipped with the Skorokhod 
topology. Then 



(S^-Sfc-.W*) :*>«}), A; = 1,2,... (4.5) 
are i.i.d. random variables in (0, 00) x D + [0, 00) with 

P(5' J7fc _ 1 - 5 ?fc _ G dy, e fc G dw) = ^e- Xy/c dyF(dw), y > 0, w G £)+[0, 00). (4.6) 
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It follows that 



{S Ck -,{e k (t):t>0}), k = 1,2, 



(4.7) 



are positioned in (— oo,a) x _D + [0,oo) as the atoms of a Poisson random measure with 

intensity c~ l \dyF(dw). Let n = max{fc > : r/ k < t q < Cfc+i}- Then S( n _ = aA9, where 
is exponentially distributed with parameter A/c. It is easy to see that 



5a-t 



for < t < S, 



for t = S, 



Cn — l 



Therefore, the first offspring in the particle system is born at time S^ n -. By (4.6) we have 
P(^Cn - S Cn - G dz) = P(e n (0) G dz) = Tl(dz), z > 0. 



By the i.i.d. property of the random variables in (4.5) we infer that {X t : t > 0} is a 
branching system with parameters (g, c _1 A, A _1 I1), where g(z) = z. In other words, the 
system have transition semigroup defined by ( 3.2[ ) and (3.8). 

Step 2) In the general case, let us consider an approximation of the subordinator with 
drift. Let {N(ds, dz)} be a Poisson random measure on (0, oo) 2 with intensity dsll(dz). 
Then a realization of {St : t > 0} is constructed by 



S t := a + 



t poo 



zN(ds, dz) — ct. 



'o 

For each k > 1 we can define another subordinator with drift {S^ : t > 0} by 

ft POO 



si 



(k) 



a + 



zN(ds, dz) — ct. 



Jl/k 



Then S 1 !^ < St and as k — > oo we have 



sup (Sl k) - S t 

0<t<T 



<?(*) 

kJrp 



S T ^0, T > 0. 



Let Tl k (dz) := l {z > 1/k} U(dz). Let {x[ k) : t > 0} be the measure-valued process defined 
by ( 4.3| with {5t : t > 0} replaced by {S^ : t > 0}. Then the first step implies that 
{Xr~t > 0} is a branching system in A^(0, oo) with transition semigroup (Qt^)t>o g iven 
by (3.11), where (u} k) ) t>o is defined as in the last step of the proof of Proposition 



3.2 



Then we can also think of {x[ k) : t > 0} as a Markov process in A^(0, oo). For t > t n > 
t n -i > ■ ■ ■ > h > and {/, /„,-■■ , fx} C C+(0, oo), we have 



Eex P {-J2(X ( t :\fd-(xi k \f)} 

i=l 

n 

= Eexp { - £<4 fc \ /,) - (X« f/S J)} 



E 



exp { - £<X«, /,) - (X| n fc) , + s k (f) 
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(4.9) 



with 



i=l 



e*(/)| < E| exp { - (Xg>, U t - t J)} " exp { - (X*\ U t %f)}\. 



(k) TT (k) 



Let (t,x) h-> TT t f(x) be the unique locally p-bounded solution of (3.6) and let (£, x) >->■ 
7r^ f(x) be the unique locally p-bounded solution of the equation with 7 replaced by 7^. 



By Proposition 3.3 and Theorem 3.4 



e k (f) <B{x£\\U t „ t J-U£lf\) 



(k) 



- U£> t J\{a) 



<*tAU t -tJ-U<j:\f\{a). 



By the proof of Proposition 



3.2 



we have U t f(x) = lim^oo U^f(x) increasingly. Then 



£k(f) — > as k — > 00. From (4.9) we get 

n 

Eexp{-^(^,/,)-(X t) /}} 

n 

= Eexp{-J2(Xt i ,fi) ~ (Xt n ,U t „ t J)}- 

i=l 

The above equality can be extended to {/, / n , • • • , fi} C B+(0, 00) by a monotone class 
argument. Then {X t : t > 0} is a Markov process in N p (0, 00) with transition semigroup 
(Qt)t>o given by (ph and (pb. □ 



5 Negative subordinators with positive drift 

In this section, we give a characterization of the branching structures in negative sub- 
ordinators with positive drift. We shall derive the result from the one in the last sec- 
tion by a time reversal approach. Suppose that II is a a-fmite measure on (0, 00) with 
J °° 1 A zll(dz) < 00. Let {S* : t > 0} be a Levy process generated by A* such that 

POO 

A*f(x)= [f(x-z)-f(x)]U(dz)+cf(x), feC\R). (5.1) 
Jo 

Assume Sq = and < c < (II, p) < 00. Then S* has Laplace exponent 

POO 

i/j(P) = c p- (l-e-^)n(dz), /3>0 
Jo 

Namely, Ee ps * = e^W. For q > let $(g) = sup{t > : ip(t) = q}. Define 

Tq := inf{t > : S? < 0}. 



14 

We have P(0 < r ~ < oo) = 1; see Corollary 5 in Section VII. 1 of pp. Then for t > 0, set 

J(t) :={ue(0,T ]:S:<t<S:_} 

and define 

ueJ(t) 

with Xg = . Note that # J(t) < oo, a.s. 

T o - 

Theorem 5.1 There is a Borel right semigroup (Qt)t>o on ^p(0, oo) defined by 

[ e-^Q t (ji,du) = e-^ u ^, /6B;(0,oo), (5.3) 

</7V p (0,oo) 

where (t, x) h-» U t f(x) is the unique locally p-bounded positive solution of 

/>t poo 

U t f(x) = /(x -*)!{.>*} + c" 1 / l {x>s} ds [l-e- u *-'^]U + (dz), (5.4) 

Jo Jo 

wiiere n + (cfe) = e _ *( ) z n(cte). Furthermore, we have 

[ (u,f)Q t (ji,du) = {ti,Tr t f), /ev3 p (0,oc), (5.5) 

JN p (0,oo) 



wiiere (t,x) (-> Tc t f(x) is the unique locally p-bounded solution of (3.6). 

Proof. Note that 0o '■= $(0) is the largest solution of V ; (/3) = 0. It follows that 

ze- zm U(dz) = V>'($(0)) > 0. 



c — 



Then (c, U + ) satisfies the conditions of Theorem 3.4 □ 



For reader's convenience, we first present a result on the distribution of time reversed 
Levy processes which should be well-known to experts. For a > 0, let {Sf : t > 0} 
and {Sf : t > 0} be two subordinators with drift starting at a > with Ee _/35 * = 
e *iW+*(°))-«0 an d Ee"^*" = e ^~ aP , respectively. Define T # (0) = inf{t > : Sf < 0} 
and T-(0) = inf{t > : 5 t ~ < 0}. Note that P{T # (0) < oo} = 1. 

Lemma 5.2 Given S*„ = a, the time reversed process {S* _ , < t < r ~} has the 
same distribution as {Sf ,0 < t < T # (0)}. 

Proof. Define I t = inf {0 A S* : < s < t} and 

J t = E 1 {5' s *</ s -}('S's - ^*-)- 
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For a > 0, set q(a) = sup{£ > : S£ — J t < x}. By Lemma 21 and Theorem 17 in Chapter 
VII of conditioned on S*_ = a, 

: < t < r } ( = } {S; - J t : < t < c(a)}. 

Then by Theorem 18 and Lemma 7 in Chapter VII of [T], under P{-|S'*_ = a}, {S* _ , : 

< t < r } has the same law as {S^ : < t < T~(0)} under P{-|T _ (0) < oo} which is 
the same as the law of {Sf : < t < T#(0)}. We have completed the proof. □ 

Theorem 5.3 The measure-valued process {X* t : t > 0} defined by Q is a single-birth 
branching system with transition semigroup {Qt)t>o determined by (5.3) and 



P{S(X*) e da} = c- 1 e-* (0)o n([a, oo))da for a > 0, 
where S(XX) = 5**_ denotes the support for XL 



(5.6) 



Proof. (5.6) follows from Theorem 17 in Section VII of [TJ. With the convention Sf_ = 
we let 

J*(t) := {ue[0,T*(0)]:S*<u<Sf}. 
For each t > define the random measure Xf on (0, oo) by 



V 



# 



E h 

ueJ#(t) 



Then by Theorem 4.1, V* is a Markov process with transition semigroup (Qt)t>o given 

we have 



by (pi). 



On the other hand, given 5** 



a, by Lemma 



5.2 



(*b -0- 



: < t < r } ( J {Sf : < t < T # (0)}. 



Thus given 5"* 



{V t * : < t < oo} = {X* : < t < oo}. 



We have completed the proof. 



□ 



6 Properties of the branching systems 

In this section we discuss the properties of the measure-valued processes via the exit 
problems for Levy processes. For a Levy process S and any x > let 

r+ = inf{t > : 3 t > x}, t~ = inf{t > : S t < x} (6.1) 

with the convention inf = oo. Set Pa;{-} = P{-|So = x }- 
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6.1 Properties of X 

In this subsection we discuss the properties of the measure-valued process X in Theo- 
rem 4.1 which is determined by process S which satisfies that Sq = a and St + ct is a sub- 
ordinator with Levy measure IT and J °° zll(dz) < c. Recall -0(A) = cA— J °°(l— e~ Xz )U(dz) 
and = sup{t > : ip{t) = q} for q > 0. Let W denote the scale function of S, i.e., 
an increasing and continuous function on [0, oo) taking values in [0, oo) with 



-Xx 



W{x)dx 



and we make the convention that W(x) = for x < 0. We will need the following solution 
to the two-sided exit problems. 



Lemma 6.1 For any t > x, y > and z > 0, 



E x e~ qT ° 



P*{r - < r t + } 



W{t - x) 
Wit) 



and 



,{S T +_ E dy, S T + - tedz, r+ < r J 
_ ^W(t-x)W(y) 



Wit) 



W(y - x))dyU(t - y + dz). 



Proof. The first identity is from the beginning of page 212 of [13]. The second identity 
follows by (8.8) of [13] with q = 0. The third identity is (8.29) of (EJ. □ 

We first present a representation of X t for any fixed t > 0. 



Proposition 6.2 The random measure X t has the same distribution as YliLo 1 w ^ 
N and (Yi) are independent random variables. 

• For a > t, 



ere 



P{N = n} 



cW(t) 



1 



cW(t) 



n > 1, 



(6.2) 



Yn = a — t and Yi, i — 1,2, 



are iid. random variables with common distribution 
ft 



l — J W(y)U(t -y + dz)dy, z > 0. 



cW(t) 

For a<t, P{N = 0} = W{t - a)/W{t) and 



(6.3) 



P{N = n} 



cW(t) 



1 



W(t - a) 



n-l 



n > 1, 



(6.4) 



(^i)i>i are z.z.a 1 . random variables with common distribution (6.3) and Y is an 
independent random variable with distribution 



W{t - a)W(y) 



W(y — a) I H(t — y + dz)dy, z > 0. 
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Proof. Observe that by the construction, the total mass X t (0,oo) is exactly the total 
number of excursions above level t, which is the same as the number of continuous down- 
crossings of level t. In addition, each excursion of S started with a jump upcrossing level 
t has to come back to level t due to overall negative drift and lack of negative jumps. 



Then (6.2) and (6.4) follow easily from the strong Markov property and Lemma 6.1 



For t < a, given N = n > 1, the excursion of S above contains n excursions at 
level t. The first excursion starts from a and all the excursions end at a. Further, by the 
strong Markov property the second to the nth excursion starts with i.i.d. initial value t + 
Yi, . . . , t+Y n -i, respectively. By the construction the support of X t is {a—t, Yi, . . . , Y n _i}. 
Note that Y\ is overshoot of the first upward jump across level t. Then by Lemma 6.1| 



P a {Vi e dz} = P t {S T + e t 
W(0) '■' 



dz, r+ < r } 



W{t) 



W(y)U(t -y + dz)dy. 



The desired result follows. The corresponding result for t > a follows similarly. 
Our next result is on the weighted occupation time for X. 

Proposition 6.3 For any f G -B+(0, oo) and h G B+(0, oo), we have 

E e -/o°° *»(*)<**./># = Ee -(^o,^o> ; 

where u is the unique nonnegative solution of the integral equation 



(6.5) 



UAX) 



-1 



l {x>s-t} 



-nds 



U(dz)[l-e 



h(s)f(x - s + t)l {x>s - t} ds. 



(6.6) 



Proof. By Theorem 5.3, similar to Section II. 3 of Le Gall [32] we can show by induction 
together with the Markov property that for any < t\ < . . . < t p and any fx, ■ ■ ■ , f p G 
5 p +(0,oo), 



where (cj t (x),t > 0, x G (0, oo)) is the unique nonnegative solution of the integral equation 

p 

H(dz)[l - e - w -«] = fi(? - u + t)M^ 



ut{x - c 



-i 



±{x>s-t}ds 



-*}■ 



Further, by taking a limit on the Riemann sums we can show that (6.3) holds. Since the 
arguments for (6.6) is similar to (5.4), one could follow the proof of Corollary 9 in Section 
II. 3 of [T7] to get (6.5). We omit the details here. □ 

It is easy to recover Laplace transform for the total occupation time (X t ,l)dt. 
Observe that it is equal to the sum of a and sizes of all the jumps of S up to time Tq , 
which is in turn equal to cTq. We then have 



E a e- q ^ {Xt ' 1)dt = E a e~ qCT o 
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6.2 Properties of X* 



Properties of the measure-valued process X* in Theorem 5.3 can also be investigated via 



the exit problems for process S* with generator (5.1 ), the negative of a subordinator with 
positive drift. 

Throughout this subsection, for q > 0, let be the scale function for the spectrally 
negative Levy process S*; i.e.; W^ q \x) = for x < and on [0, oo), it is an increasing 
and continuous function taking values in [0, oo) with 

OO -I 

e- Xx W (q \x)dx 



V(A)-g' 



for A > := sup{A > : if)(\) = q}, where i[)(X) = cA - / °°(1 - e- Xz )U(dz). Write W 
for W^°K We also first present a result on the two-sided exit problems of S*; see Chapter 
8 of [13] and 0. 

Lemma 6.4 For any y > > z, 

Ee~ qT ° = 1 - 4r^ (5) (0) 
$(g) 

and 

P{S*-_ E dy, S*- G dz) = W(0)e- my U(dz - y)dy. 

Proof. The first identity is just (8.6) of [13J. The second identity the equation right 
after (8.29) in [Lj. □ 

For any y > t, let 

W{z) 



and 



g(y) := c -^^ J U(y - dz) 
h(y) := c-^TOM | jT H(y - (l - + n((y, oo)) 



One will see from the proof of Proposition 6.5 that J t °° g(y)dy + J f °° h(y)dy = 1. 

Fix t > until the end of the following Proposition |6.5| We first proceed to recover 
distribution for the total mass for . The proof of the following representation result is 



similar to Proposition 6.2 and is omitted 



Proposition 6.5 The random measure X\ has the same distribution as ^2f =Q l where 
N and (Yi) are independent random variables. 

' ~ cW(t) 



P{N = 0} = 1 - — \— (6.7) 



and for any n > 1 

i / poo \ n—1 r-OO 

P{N = n} = cWttj{J t 9{y)dy J I Hy)dy - (6i 
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Y has the density function h(t + y)/ L h{r)dr, y > and Yi, i = 1, 2, . . share the 
common density function g(t + y)j f g(r)dr, y > 0. 



Proof. Since N = if and only if the whole excursion of S* stays below level t up 
to time t ~ , the probability (6.7) just follows from Lemma 6.1 Observe that the total 
mass X^(0, oo) is exactly the number of up-crossings (the same as the number of down- 
crossings) of level t by process S* until the time Tq . Each up-and-down-crossing of level 
t corresponds to an excursion starting at level t. All of such excursions end at level t 
except that the last one ends below at time t ~ , where the last excursion determines 
the residual life time of a particle that can be either the ancestor or an offspring. Using 
solutions to the two-sided exit problem in Lemmas 6.1 and 6.4 together with the strong 
Markov property repeatedly at those up-crossing times of level t we have 

P{X*(0,oo)=n} 

= P{r+ < r -} (7V{S;- G dz}P z {r+ < r -}) 

x (J ?t{S* rr e dz}P z {r+ > r } + P t {S* r - < 0}) 



W(0) 
W{t) 

x 



n-l 



c -»(o)(^*)W(0)dy|jf n(y-dz)(V^)+n((y,oo))}. (6.9) 

Therefore, the probability (6.8) follows. 

Given X*(0, oo) = n, the support of X*(0, oo) consists of those distances between the 
pre-down-crossing (of level t) values of S* and t for the n excursions from t. By the strong 
Markov property all these distances are independent. By Lemma 6A the distances for the 
first n — 1 excursions following the same distribution of 

/ P t {S* r -_ et + dy, S* r - G dz}P z {r t + <r }( [ P t {S* T - G dz}P z {r+ < r } 
Jo 1 ' \Jo * 



e -*WvW{0)dy / U(t + y-dz) 



W{z) 



g(r)dr 



g(t + y)dy 



g(r)dr 



The distance for the last excursion follows the distribution of 



Qf Pt{S* r -_ et + dy, S* r „ G dz}P z {r t + >r } + P t {S* r -_ et + dy, £T_ < 0}^) 
x £ P t {S* T -_ et + dy, S* r - e dz}P z {r+ > r } + P t {S* r -_ et + dy, 5*_ < 0}^) 

= h(t + y)dy I / h(r)dr 



-i 
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Our next result is on the weighted occupation time for X*. The proof is similar to 



Proposition |6.3| and is omitted. 

Proposition 6.6 For any f G S+(0, oo) and h G B+(0, oo), we have 

E e - h (t)(x?J)dt _ g e -(x *. w o) ) (6.10) 

where u is the unique nonnegative solution of the integral equation 

poo poo 

co t (x) - C l / l {M} ds / Yi + {dz)[l - e- w 'W] 

</i JO 

h{s)f{x - s + t)l {x>s _ t} ds. 



Observe that the total occupation time (X^ , l)dt is just the sum of the sizes of all 
the jumps of process S* before time r Q ~ together with S*-_. Further, this sum is equal to 



ct since Sq = 0. By Lemma 6.4 we then have 



Eg-iJ-^v)* = Ee~ qCT ° = 1 ^-W {qc) (0) = 1 ^— . (6.11) 

$(gc) $(gc) 

7 Connections with the CMJ model 

Informally, the Crump-Mode-Jagers branching processes or the CMJ process counts the 
size of a branching population system with random characteristics. Informally, a particle, 
say x, of this process is characterized by there random process 

Cc(')j u x) 

which is an i.i.d. copy of (A, £(•), oo) and the reproduction scheme is given in the following 
sense: if x was born at time a x , then 

1. X x is the life length of x\ 

2. ( x (-) = {0 < ( x < ( x < ■ ■ ■ < \ x } is a point process defined on (0, X x ). {Q + o~ x : 
i — 1, •} is the collection of splitting times of x at which it produces offspring. 

3. u) x is the number of children produced by x at time cx^ + ( x - 

Let denote the total number of individuals in the system at time t with Z(0) ances- 
tors. In general, the process {Z(t) : t > 0} is not Markovian unless X x is exponentially 
distributed. Now assume that 

1. The distribution of A is determined by a probability measure r](dx) on (0, oo); 

2. £(•) is a Poisson point process with parameter a; 

3. The distribution of uj 1 is determined by a generating function g(-). 
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According to the argument in Section 2 and j3], we may define a measure- valued Markov 
process Y = {Y(t) : t > 0} with transition probabilities given by 

/ e-^Q t {fi,diy) = e-^\ /e5 + (0,oo), (7.1) 

JN(0,oo) 

where (t,x) i-» U t f(x) is the unique locally bounded positive solution of 

U t f(x) = f(x-t)l {x>t} + a [ l {x>t „ s} [l-g(( V ,e- u *f))]ds. (7.2) 

Jo 

Then the CMJ process {Z(t) : t > 0} is just the total mass process of Y; i.e. Z(t) = 
(Y(t),l). 

The connection between Levy processes and CMJ processes was first investigated by 
Lambert in [14] which showed that the contour process of a splitting tree defined from a 
suitable CMJ process is a spectrally positive Levy process with negative drift killed when 
it hits 0. The starting position of the Levy process is just the life time of the ancestor. 
Equivalently, given such a Levy process, one could construct a CMJ process; see also [T3] . 
In those works, the Levy measure, say 7, is assumed to be a cr-finite measure on (0, 00] 
with J( 0oo ] 1 A z r y{dz) < 00. Our main result, Theorem 3.2, also gives similar relationships 



between one-sided Levy processes of bounded variation and CMJ processes. 
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